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A single equation governing the various properties encountered in propa- 
gation of Gaussian beams is discussed. These characteristics may be graph- 
ically presented on an impedance chart such as a Smith chart or in the form 
of normalized curves. The geometrical representations highlight the relation 
between Gaussian mode theory and geometrical optics formulas. 

The power coupling coefficient between two Gaussian modes suggests a 
complex mismatch coefficient whose geometrical representation is essentially 
the same as that of a complex reflection coefficient in transmission-line theory. 
Application of the Smith chart in determining a complex mismatch coeffi- 
cient is illustrated by graphical solution of a beam-matching problem. 

I. INTRODUCTION 

The propagation of a Gaussian beam and its transformation through 
a lens has been well treated in previous literature. 1 •*•'•* This paper will 
show that a single, formally identical equation governs three properties 
of Gaussian beam propagation: (i) the phase front curvature and the 
beam radius in terms of the distance from the beam waist and the 
minimum beam radius; («) the propagation of a Gaussian beam in free 
space; (Hi) the transformation of a Gaussian beam through a lens. 
Geometrical representations of these characteristics highlight the re- 
lationship between Gaussian beam propagation and geometrical optics. 

Several recent papers 6,6,7,8 have been devoted to graphical solutions 
of Gaussian mode problems. One recalls that the Smith chart is a geo- 
metrical representation of complex reflection coefficient in transmission 
line theory. It seems logical, therefore, to look for the counterpart of a 
complex reflection coefficient in Gaussian mode theory in order that the 
full potential of the Smith chart may be realized in graphical solutions 
of Gaussian mode problems. In this paper, a complex mismatch coef- 
ficient will be defined such that the geometrical representation of this 
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coefficient is essentially the same as that of a complex reflection coef- 
ficient in transmission line theory. This observation immediately sug- 
gests the application of the Smith chart of complex mismatch coefficients 
to a graphical solution of beam matching problems. 

II. THE ANALOGY AMONG THREE PROPERTIES OF A GAUSSIAN BEAM 

In order to discuss the Gaussian beam transformation it is convenient 
to summarize the relationships among the parameters of a Gaussian 
beam propagating along the z-axis first. 4 ' 9 

w(2) = • 1/ 1 + (£)" (la) 

' w -[ l + (£)■]■ (lb) 

In (1), w is the beam radius at which the field amplitude has fallen to 
1/e of its maximum value on the z-axis, w is the minimum beam radius 
(called the beam waist) where one has a plane phase front at z = 0, and 
r is the radius of curvature of the phase front at z. It should be noted 
that the phase front is not exactly spherical; therefore, its radius of 
curvature is exactly equal to r only on the z-axis. The z-coordinate, 
measured from the beam waist, is taken to be positive to the right and 
negative to the left; the parameters are illustrated in Fig. 1. Equation 
(1) may be solved for w and z as follows: 

_ w 

■■■*- (2a) 
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2 ~ 1+ /^Y '-"■■ 






Equations (1) and (2) may be transformed into a single equation of 
complex variables: 

1 -,-i. ' , „ . (3) 



irw 2 /\ r (tw 2 /\) + iz' 

We are also interested in the transformation of a Gaussian beam going 
through a lens of focal length/ as shown in Fig. 2. A beam with its mini- 
mum beam radius Wi located at rfi will become another beam with its 
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Fig. 1 — Parameters of a Gaussian beam. 

minimum beam radius iv 2 located at (h . Since the beam radius remains 
the same in passing through the lens we obtain from (la) 



Wi Vl + (Xdi/wilh 2 ) 2 = w> 2 Vl + (M./widi 1 ) 2 . 



(4) 



The thin lens formula states that the change of the phase front curva- 
ture may be approximated by the reciprocal of the focal length. Thus, 
using (lb), one obtains 
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Fig. 2 — The analogy between the Gaussian beam transformation through a lens 
and the Gaussian beam propagation in free space. 
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+ 77; i / -2 AJ x2, = 7 (5) 



di[l + (xW/Adi) 2 ] ' d,[l + (TtfcVXdtn / 

where the sign of the second term on the left side is positive because d 2 
has been taken to be positive to the left of the beam waist following the 
convention of Goubau 1 and Kogelnik. 2 Straightforward algebra will lead 
to the following solutions for (4) and (5) : 

d * - i = t (6) 

J (di , \ , (irWi 



( f -,j + c 



\f 



* _ iY + ^"v (7) 



The above two expressions are essentially rearrangements of formulas 
obtained by Goubau. 1 When iriui 2 /\ « | di — f \, (6) and (7) approach 
the thin lens formula and the magnification formula of geometrical 
optics. * When TWi 2 /\f — > °° the condition for geometrical optics focusing 
of parallel rays is obtained. These two equations have been plotted in 
Figs. 3 and 4 for various values of p = irwi 2 /\f. There it is seen that the 
singularity at di = / in geometrical optics is eliminated in Gaussian mode 
transformations. The maxima and minima in Fig. 3 may be easily found 
by differentiating (6) with respect to di/f, and they are given by 

d 2 , . / i di irwi / 8 \ 

- = 1 ± _ :... when -j = 1 ± -— r . (8) 

/ 2tWi 2 /\ f A/ 

The points of inflection are 

dl =l and liV^' 



Equations (6) and (7) can be combined into one equation of complex 



variables : 



(f)"< 



- i 1 1 2 - 1 1 = 



7 r (t>ih>) (9) 

* Simple algebra may easily reduce them to the more familiar forms, 3" + T = 

* and -z-. = j-r . 

' Wi 2 di 2 
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Fig. 3 — The graphical presentation of (5) and (14). 



Next, let us consider the propagation of a Gaussian beam along the 
2-axis. One would like to express the beam radius w 2 and the curvature 
l/j'2 of the phase front at any given point in terms of the beam radius w\ 
and the curvature \/i\ at some other point along the axis. Since any 
Gaussian beam has a single beam waist at a definite location, one obtains 
from (2) the following conditions: 



uh 



W-2 



Vi + (WAn) 2 

r,. 



y/l + (7rw 2 2 /Xr 2 ) 2 



n 



i i (\ I 2\2 — #2 — Z\. 

1 + (Xn/TWi ) 



1 + (Xr2/irWo ) 
Let L = zi — Z\ . Rearranging (10) and (11) yields 



(10) 
(ID 
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irWi /X ttw-i /X 



wwi /X/ \n/ \tw 2 /X 



+ " 



© 



(12a) 



and 



1_ 1_ 



i v + l, , - l + 



If we also rearrange (4) and (5) as follows : 

7riyi 2 /X 7rti)2 2 /X 



a 



= L. (12b) 



(xWVx) 2 + di 2 " (W/x) 2 + di 



(13a) 



"a + / - 2/. \2 , ,2 = 7 (lob; 



(WA) + di (iriJk/\) + d 2 / 

the analogy* between (12) and (13) now emerges with the following 
one-to-one correspondences : 
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Using the above analogy, one may immediately write down the solu- 
tions of (12) for w> 2 and l/r 2 

k - i = d (H) 

I _ L _ -A . I L 



r\ / \jrWi 2 /\ 



* The analogy would look even better, if n had been assumed positive to the 
left of the beam waist. This fact, of course, results from the sign convention of di , 
which is positive to the left of the lens. 
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L 



L/ 



irw\ 



irw-r/X 






+ 



L 



irtVi 2 /\ 



(15) 



Except for slightly different normalizations, these two expressions coin- 
cide with two formulas obtained by Rowe 10 using rather involved alge- 
bra. The curves in Figs. 3 and 4 represent (14) and (15), if the coordi- 
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Fig. 4 — The graphical presentation of (6) and (15). 
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nates and the parameters are replaced by those shown in parentheses in 
the figures. When the wavelength approaches zero, the Gaussian beam 
becomes the spherical wave of geometrical optics. The singularity of the 
spherical wave is eliminated in the Gaussian beam representation. 
Equations (14) and (15) can also be combined into one equation of 
complex variables: 

-<i'5-iU 



wtf/X \r 2 ) ( L \, i (_L_{\ (16) 



It is interesting to observe that (3), (9), and (16) are formally identical, 
and they are amenable to graphical representation on an impedance 
chart, e.g., a Smith chart. Thus, the circumference of the unit circle on 
a Smith chart corresponds to geometrical optics, while the interior de- 
scribes the properties of a Gaussian beam. The transformation may be 
performed by simply taking diametrically opposite points on a Smith 
chart. It should be noted that while (3) has not been normalized, (9) 
and (16) are dimensionless. Representation of (3) by the Smith chart has 
also been proposed by Deschamps and Mast; 7 the analogy between (9) 
and (16) corresponds to the dual forms of the cartesian Gaussian beam 
chart. 66 Rearrangement of (14) and (15) will yield equations similar to 
those discussed in Ref. 8. 

III. THE COMPLEX MISMATCH COEFFICIENT DIAGRAM 

Using some of the results of the previous section, a graphical solu- 
tion of mode matching problems using a properly defined complex mis- 
match coefficient will now be discussed. Equation (3) can be rewritten as 

..2 J 



irW 

— IZ = 



W 2 , 



r ~ / i \ , .i d9) 

+ 1- 



- 2 /a/ ' r 
and one may adopt the following notations : 



irw 2 /\ 



+ i- = R + iX = Z (20a) 



^ -iz = G -iB = Y. (20b) 

A 

Thus, a Gaussian beam may be characterized by either a beam impedance 
Z or a beam admittance Y. Now the following interesting identity can 
easily be verified. 
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Z\ — Zo 



Zx + Z„* 



Y x - F 



|2 



= 1 - T (21) 



where the subscript 1 represents the incoming beam, represents the 
fundamental mode of a receiving system,* denotes complex conjugate, 
and t is the power coupling coefficient 11 between two Gaussian modes. 
The first two expressions in (21) certainly look like reflection coefficients 
in transmission lines. This observation suggests the designation of 
(Zi - Z )/(Zi + Z„*) and (Fi - F )/(Fi + F *) as complex mis- 
match coefficients. They are indeed identical to reflection coefficients in 
form except for the complex conjugate in the denominator. If Z\ (or Fi) 
is normalized with respect to a real Z (or a real F ), then the complex 
mismatch coefficient is formally identical to a complex reflection coef- 
ficient. Both the complex mismatch coefficient and the complex reflection 
coefficient are bilinear transformations of the impedance or admittance. 
The Smith chart is a geometrical representation of the complex reflec- 
tion coefficient. The same Smith chart may also represent the complex 
mismatch coefficient. With respect to a reference mode, any phase front 
of a Gaussian beam may be represented by two diametrically opposite 
points on a Smith chart, i.e., the normalized admittance and the nor- 
malized impedance. The propagation of a Gaussian beam corresponds to 
travel along a conductance circle, while imposing change of radius of 
curvature by a lens, say, corresponds to travel along a resistance circle. 
The pair of two diametrically opposite points may be reduced to one 
point by imposing a flip-over Smith chart on the original Smith chart, 
however, overcrowded coordinate lines are not desirable in practical 
graphical solutions. This latter version of geometrical representations 
corresponds to a bilinear transformation of the Gaussian beam chart 
discussed in Ref. 5, and is the same as that suggested in Ref. 7. It 
should be noticed that Collins' Gaussian beam chart 5 is a geometrical 
representation of the beam impedance, while Li's dual form 6 is that of 
the beam admittance. Here, however, the Smith chart is utilized as a 
complex mismatch coefficient diagram and its relationship to the power 
coupling coefficient between two Gaussian modes is identified. 

IV. A NUMERICAL EXAMPLE OF BEAM MATCHING 

In order to illustrate the application of the Smith chart, consider the 
problem of matching the output beam of an optical maser to an inter- 
ferometer as shown in Fig. 5. The maser resonator consists of mirrors of 
5 and 10-meter radius of curvature separated by 1 meter, and the inter- 
ferometer consists of a pair of 10-meter mirrors separated by 10 cm. 
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Fig. 5 — The geometry of a beam matching problem. 

The maser beam exits from the 5-meter mirror and a lens with a focal 
length of 2 meters is available for matching. First, one calculates the 
minimum beam half-widths and their locations for both resonators by 
the following formulas: 4 



(22) 



■kw Vd(Ri - d)(R 2 - d)(Ri + fl 2 - 


- 2d) 


X Ri + R 2 - 2d 

„ _ , Rz — d 

Ri -f- Rz — 2d 





(23) 



where R\ and Ri are the radii of curvature of the mirrors, d is mirror 
separation, w is the minimum beam half-width, and S is the location of 
the beam waist from the mirror of radius Ri . Now we readily get 
7rtDi 2 /X = 0.705 m, Si = 0.05 m for the interferometer, and irw£/\ = 
1.727 m, $2 = 0.692 m for the maser resonator. If all the parameters are 
normalized with respect to tWi 2 /\, then the fundamental mode of the 
interferometer is represented by the unity conductance circle passing 
through the center of Smith chart in Fig. 6 where corresponds to the 
beam wasit for the interferometer. The maser output beam is represented 
by the 2.45 (= w-f/wi 2 ) conductance circle passing through M which 
corresponds to the beam waist for the maser resonator. The inversions of 
these two circles with respect to the center give the beam radius and 
the phase front curvature along the propagation path of the beam. 

Let it be required that insertion of a 2-meter focal length lens changes 
the phase-front curvature of the maser output beam to that of the in- 
terferometer beam at the point where the two beams meet. Since the 
beam radius remains constant going through the lens, application of the 
thin lens formula demands that a segment of the beam resistance circle 
between the inverted maser and interferometer circles be equal to the 
length A^B' = rio^/Xf = (rw^/X^l/n - l/r 2 ) = 0.3523. After deter- 
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Fig. 6 — The graphical solution of a beam matching problem. 



mining the two pairs of diametrically opposite points A-A' and B-B' 
(Fig. 6), the distances from the lens to the interferometer beam waist 
and to the maser beam waist can be read off the chart as OB = 4.3 X 
0.705 = 3.03 m and MA = 6.5 X 0.705 = 4.58 m. The distances fiom 
the lens to the interferometer input mirror and to the maser output 
mirror are 3.03 - 0.05(3,) = 2.98 m and 4.58 - 0.69(S 2 ) = 3.89 m, 
respectively. 
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V. DISCUSSION 

The advantage of this complex mismatch coefficient diagram over the 
Gaussian beam charts discussed in Refs. 5 and 6, is similar to that of the 
Smith chart over cartesian impedance and admittance charts. Neither 
of the latter can represent an arbitrary transformation of aGaussian beam, 
whereas the present mismatch coefficient chart, at least in principle, 
encompasses all possible transformations. Curvilinear coordinates are 




Fig. 7 — An example of approximate beam matching. 
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not a problem here because a good network of circles is readily available 
in the form of Smith charts. Sometimes it may be impossible to achieve 
perfect beam matching due to limitations imposed by the inadequacy of 
a lens or of space. One is able to minimize the mismatch using the 
procedure discussed here. 

In the above example, there exists a minimum focal length of the lens 
beyond which perfect matching is not possible. In Fig. 6, this minimum 
is represented by the maximum arc length CD' (0.64) of the beam 
resistance circle between the inverted maser and interferometer circles 
which corresponds to a focal length of vwf/XCD' = 0.705/0.64 = 1.102 
m. If one had only a one meter focal length lens available, CD' would 
be extended to E' as shown in Fig. 7 and the reciprocals of C and E' 
are then found to be C and E. The distances from the lens to the maser 
output mirror and to the interferometer input mirror would be CM — 
Si - 1.8 X 0.705 - 0.05 = 1.22 m and EF - S 2 = 1.56 X 0.705 - 
0.69 = 0.41 m for a resulting mismatch loss \ OF \ 2 = 0.015. 
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